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Abstract 

Factorizable amplitudes in B decays to light vector meson pairs give a longitudinal 
polarization satisfying 1 — /l = 0{l/rn^). This remains formally true when non- 
factorizable graphs are included in QCD factorization, and is numerically realized 
in B — > pp. In AS = 1 decays a QCD penguin annihilation graph can effec- 
tively contribute at leading power to the transverse and longitudinal amplitudes. 
The observed longitudinal polarization, fi{B — * 4>K*) ~ 50%, can therefore be ac- 
counted for in the SM. The ratio of transverse rates r_i_/T|| provides a sensitive test 
for new right-handed currents. The transverse b — > sg dipole operator amplitudes 
are highly suppressed. CP violation measurements can therefore discriminate be- 
tween new contributions to the dipole and four quark operators. SU(3)p violation 
in QCD penguin amplitudes can easily be 0(1), in general, due to annihilation. 
Implications for B — > pK* polarization and New Physics searches are pointed out. 
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1 Introduction: 'helicity-flip' suppression 



Polarization in B — > VV decays should be sensitive to the V— A structure of the Standard 
Model due to the power suppression associated with the 'helicity-flip' of a collinear quark. 
For example, in the Standard Model the factorizable graphs for B — ► <pK* are due to 
transition operators with chirality structures (sb)v-A(ss)v^A, see Figure. 1. There are 
three helicity amplitudes, A , A~, and A + , in which both vectors are longitudinally, 
negatively, and positively polarized, respectively. In A~ a collinear s or s quark with 
positive helicity ends up in the negatively polarized <f), whereas in A + a second quark 
'helicity-flip' is required in the form factor transition. In the case of new right-handed 
currents, e.g., (sb)v+A(ss)v±A, the hierarchy is inverted, with A + and A~ requiring one 
and two 'helicity-flips', respectively. 

Helicity-flip suppression can be estimated by recalling that the probability for a 
positive helicity free fermion to have negative spin along some axis is given by sin 2 0/2, 
where 9 is the angle between the axis and the momentum vector. For a <fi meson in a 
symmetric configuration the transverse momentum of the valence quarks is k± ~ m^/2, 
implying that the helicity suppression in A- is ~ m^/ms. The form factor helicity 
suppression in A+ should be approximately pr/mb, where pr is the transverse momentum 
of the outgoing s quark. The latter can be estimated by identifying it with the transverse 
momentum of the b quark. In the 'Fermi momentum' model of P < p\ > x I 2 k, p F j^?>. 
Using the equivalence of this model to a particular HQET based shape function ansatz 
and for illustration taking A « 500 MeV and -Ai » 0.3 GeV 2 yields p F » 400 MeV, 
or a helicity suppression of ~ 0.05. 

These simple estimates should be compared to naive factorization, supplemented by 
the large energy form factor relations [3] (also see jl]). For B — > <pK*, 

A = i^X s t ~aA%, A^ = i^\l~aAl^. (1) 

The coefficient a = a 3 + a 4 + a 5 — |(a 7 + a 9 + a 10 ), where the a« are the usual naive 
factorization coefficients, see e.g. [3], and = V p bV* q . The large energy relations imply 

A VrV 2 = fv 2 m 2 B Cl\ A v^v 2 = -fv 2 mv 2 mB2Cl 1 , A$- lVa = -f V2 m V2 m B 2(l 1 r l [ 1 . (2) 

We use the sign convention (V|g7 M g|0) = — i/ym^f*. £|| and (± are the B — > V form 

factors in the large energy limit Both scale as m b 3 ^ 2 in the heavy quark limit, 
implying that helicity suppression in A~ is ~ m^/ms which is consistent with our 
estimate (the form factor transition contributes 2(± in Ay iV2 ). r±_ parametrizes the form 
factor helicity suppression. It is given by 

_ (1 + m v Jm B )A v l 1 - (1 - m Vl /m B )V v ^ 
(l + m Vl /m B )AX x + (1 - m v Jm B )V v ^ ' 

where A\^ and V are the axial-vector and vector current form factors, respectively. The 
large energy relations imply that it vanishes at leading power, because helicity suppres- 
sion is 0(l/m). Light-cone QCD sum rules [BJ, and lattice form factor determinations 
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Figure 1: Quark helicities (short arrows) for the Standard Model naive factorization 
coefficients a^. Upward lines form the emission vector Vi- q = s for B — > <pK* . 

scaled to low q 2 using the sum rule approach [7j, give r*[* rs 1 — 3 %; QCD sum rules give 
r*[* ~ 5% 8J; and the BSW model gives rf* « 10% jH]. These results are consistent 
with our simple estimate for form factor helicity suppression. 

The large energy relations giving rise to (J2J) are strictly valid for the soft parts of the 
form factors, at leading power and at leading order in a s . However, the soft form factors 
are not significantly Sudakov suppressed in the Soft Collinear Effective Theory (SCET) 
[TU] . The results of jU[TT] thus imply that the 0(a s ) form factor contributions, particu- 
larly the symmetry breaking corrections to the large energy relations, can be neglected. 
In fact, r± does not receive any perturbative corrections at leading power [121 IH ITT] : 
again, this is because form factor helicity suppression is 0(l/m). Furthermore, power 
corrections to all of the form factor relations begin at 0(1/ m) (rather than 1/y/m) in 
SCET ^3]. Therefore, the above discussion of helicity suppression in naive factoriza- 
tion will not be significantly modified by perturbative and power corrections to the form 
factors. 

In the transversity basis [TJ] the transvese amplitudes are A±,\\ = (A~ =F A + )/V2 
(A±,\\ = (A + =F A~)/y/2 ) for B (B) decays. The polarization fractions satisfy 

1 - f L = O {l/ml) , /j.//|, = l + 0(l/m b ), (4) 

in naive factorization, where the subscript L refers to longitudinal polarization, /j = 
Fi/r total , and fL + f± + f\\ = 1- The measured longitudinal fractions for B — » pp are 
close to 1 ^SIQSIEI- This is clearly not the case for B — * <f)K*°, for which full angular 
analyses yield 

fi, = .43 i .09 i .04, f ± = .41 ± .10 ± .04 [THj (5) 
fi, = .52 i .07 i .02, f ± = .27 ± .07 ± .02 [JJj|. (6) 

Naively averaging the Belle and BaBar measurements (without taking large correlations 
into account) yields f±/f\\ = 1.4 ± .7. In the charged mode, BaBar has measured 
fi((j)K* + ) = 0.46 ± 0.12 ± 0.03 JE]- We must go beyond naive factorization in order to 
determine if the small values of fi((f>K*) could be due to the dominance of QCD penguin 
operators in AS* = 1 decays. 

2 QCD factorization for B — > VV decays 

In QCD factorization [20] exclusive two-body decay amplitudes are given in terms of 
convolutions of hard scattering kernels with meson light-cone distribution amplitudes. 
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At leading power this leads to factorization of short and long-distance physics. This 
factorization breaks down at sub-leading powers with the appearance of logarithmic in- 
frared divergences. Nevertheless, the power-counting for all amplitudes can be obtained. 
The extent to which it holds numerically can be checked by introducing an infrared 
hadronic scale cutoff, and assigning large uncertainties. Non-perturbative quantities are 
thus roughly estimated via single gluon exchange. In general, large uncertainties should 
be expected for polarization predictions, given that the transverse amplitudes begin at 
0(l/m). However, we will find that this is not the case for certain polarization observ- 
ables, particularly after experimental constraints, e.g., total rate or total transverse rate, 
are imposed. Our results differ substantially from previous studies of B — > VV in QCD 
factorization [231 El] ■ Of particular note is the inclusion of annihilation topologies. The 
complete expressions for the helicity amplitudes are lengthy and will be given in [22]. 
Expressions for a few contributions are included below. 

In QCD factorization, the Standard Model effective Hamiltonian matrix elements can 
be written as [2T1I22] 

{v^n^^Y.K^M^+rs^B) , (7) 

" p=u,c 

where h labels the vector meson helicity, and D = s(d) for AS = 1(0). Tg gives rise to 
annihilation topoplogy amplitudes, to be discussed shortly, and 

10 
i=l 

where p — u,c, and V[V{ = V1V2 or V2V1. The coefficients a^ ,h contain contributions 
from naive factorization, vertex corrections, penguin contractions, and hard specta- 
tor interactions. The transition operators j{ <g> j\ are 5 pu (ub)v-A <8> (Du)y-A 0=1) j 
5 pu {Db)v- A ® (uu)v- A (i=2); [{Db) V -A ® (qq)v T A (i=3,5)] xfe, (i=9,7); \{qb)v-A ® 
{Dq) v _A (i=4)] x|e g (i=10); [(qb) S -p ® {Dq) s +p (i=6)] x|e g (i=8), where q is summed 
over u,d,s. For i=l,2; 3-6; and 7-10 they originate from the current-current Qi,2] QCD 
penguin <53,..,6j an d electroweak penguin operators Qr,..,io, respectively. For i ^ 6, 8, 
(V{ h V2 h \j{ ® j|p) is defined at leading-order as 

(V! h \j{\B)(V 2 %m = -ic\V{Vi) A h v , v , (9) 

The c coefficients contain factors of ±1, ±l/y2, arising from the vector meson flavor 
structures. V 2 ' (V[) is the 'emission' ('form factor') vector meson, see Figure 1. The 
i = 6, 8 matrix elements vanish at tree-level, i.e., at leading order in a s , as local scalar 
current vacuum-to- vector matrix elements vanish. Due to the underlying flavor structure, 
the effects of a 3 -a 10 are describable in terms of a reduced set of coefficients [22] 

h h 1 h p,h p,h . fp,h flfA 

a 3(3EW) — a 3(9) + a 5(7)> a 4(4EW) — a 4(10) + a 6(8) ' \ W ) 
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where dg^ - ou, 6 ^ 
matrix elements in a 



and (j±' 8 <S> j^} are next-to- leading order 



in which j 2 again forms the emission particle V%. The arguments 
(V1V2) are understood throughout. 



At next-to- leading order, the coefficients a[ h can be written as [22 



(W 



a 



4vr 2 



V?(K 2 ) +—H?{VM 



+ P t h ' p (V 2 ), (11) 



where the upper (lower) signs apply when i is odd (even). The superscript 'p' appears for 
% = 4,6,8, 10. The iVj are tree-level naive factorization coefficients (Figure 1); at next- 
to-leading order the V/ 1 account for one-loop vertex corrections, the Pf' h for penguin 
contractions (Figure 2), and the for hard spectator interactions (Figure 3). They 
are given in terms of convolutions of hard scattering kernels with vector meson and B 
meson light-cone distribution amplitudes. For each i, the corresponding graphs have the 
same quark helicity structure. 

Two twist-2 light-cone distribution amplitudes <f>\\(u) and (p±(u), and four two-particle 
twist-3 distributions (and their derivatives) enter the longitudinal and tranverse vector 
meson projections j2E]- The argument u (u = 1 — u) is the quark (antiquark) light-cone 
momentum fraction. The two-particle twist-3 distributions can be expressed in terms 
of 0||,_l(w) via Wandura-Wilzcek type equations of motion if higher Fock states are 
ignored. The twist-3 vector meson projections then depend on the three distributions, 



, u 



= I dv 



= I dv 
'0 



Hi* 



; u 



= I dv 
'0 



<j)j_{v) 



dv 



(12) 



$ a and $b project onto transversely polarized vectors in which the quark and antiquark 
flips helicity, respectively. $„(«), defined in [22] , projects onto longitudinally polarized 
vectors in which either the quark or antiquark flips helicity. Light quark mass effects 
are ignored, and a discussion of twist-4 distribution amplitudes and higher Fock state 
effects is deferred j2H]- The leading-twist distribution amplitudes are given in terms of 
an expansion in Gegenbauer polynomials [2*61 IE]. 



4>i(u, jj) = 6uu 



n=l 



i =±, 



(13) 



Our numerical results include the first two moments ati a 2} i- The asymptotic forms of 
the twist-3 distribution amplitudes are 3m 2 , 3m 2 , and 3(w— u) for $ (m), $b(w), and $ v (u), 
respectively. The B light-cone distribution amplitude 0+ jl] enters the hard spectator 
interactions through its inverse moment, j dlcj)^/ I = itib/^b- 

The naive-factorization coefficients are iVj = (i=6,8), iV, = 1 (i^6,8 ). The longi- 
tudinal (transverse) helicity amplitudes arise at twist-2 (twist-3) in naive facotrization. 
The vertex corrections are negligible compared to the theoretical uncertainties. Note 
that Vgg 7^ for all h. At 0(a s , a) the longitudinal penguin contractions are, respec- 
tively, P 4 °;fo(^) = P*,w{V) at twist-2 and P^(V) = - P 6>8 (V) 2 m v J^)/m h ^)f v at 
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Figure 2: Quark helicities for Ff' : Qi^^s (left), Q% g (right). For P§ ,h flip the q,q 
helicites. Upward lines form V2. 

twist-3. The quantities Pi(V) are the VP counterparts defined in [22j, and fy is the 
scale-dependent tensor-current decay constant. The transverse penguin contractions are 
P§'% = to twist-4, and at twist-3, 



P^ P (V 2 ) = (a s C F /4 7 riV c )(C 1 [ J L-G±( Sp )]+C3[2L-G|(0)-G'±(l)] 



± 



+ (C 4 +C 6 )[L-3G±(0)-G±( Sc )-G±(l)]), 
= (a/9n) (d/JV c + C 2 )[L - G± ( Sp )], 



(14) 



where L = (4/3) lnm^/ju + 2/3, L = (20/3) lnmj/ju, s p = rn 2 v jrn\ (s u = 0, s c = m 2 c /ml 



du G(s — ie,l — u) §1\{ U ) , 



(15) 



and G(s,x) is the well known penguin function, see e.g., [21j . The penguin contractions 



account for approximately 30% and 20% of the magnitudes of a% and a c A ~ (for default 
input parameters), respectively, before including the hard spectator interactions. 

The dipole operators Q$ g , Q77 do not contribute to the transverse penguin amplitudes 
at 0(a s ) due to angular momentum conservation: the dipole tensor current couples to a 
transverse gluon, but a 'helicity-flip' for q or q in Figure 2 would require a longitudinal 
gluon coupling. Formally, this result follows from the Wandura-Wilczek relations and 
the large energy relations between the tensor-current and vector-current form factors 
|25j . For example, the integrand of the convolution integral for P 4 ~ vanishes identically, 
whereas P^ oc Jdu(^ 2 (u) — $^ 2 (m)) = 0. Note that transverse amplitudes in which a 
vector meson contains a collinear higher Fock state gluon also vanish at 0(a s ) . This can 
be seen from the vanishing of the corresponding partonic dipole operator graphs in the 
same momentum configurations. Transverse 0(a 2 s ) spectator interaction contributions 
are highly suppressed and are studied in |23j . 

The hard spectator interaction quantities H\ contain logarithmically divergent in- 
tegrals beyond twist-2, corresponding to the soft spectator limit in Vj, see e.g., Figure 
3. We integrate the quark light-cone momentum fraction in V\ over the range [0, 1 — e], 
and replace the divergent quantities \ne with complex parameters Xh. As in |2*H I22j. 
these are modeled as X H = (1 + g H e lipH )\rLm B /A h , with qh < 1 and A h f& 0.5 GeV. 
This reflects the physical (9(Aqcd) infrared cutoff, and allows for large strong phases 
from soft rescattering. For z^6,8: H® first arises at twist-2; arises via a twist-3 V2 x 
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twist-2 Vl projection; and Hf arises via a (twist-3) 2 projection. For i = 6, 8: Hq 8 = to 
twist-3; Hq 8 arises via a twist-3^ 1 x twist-2 V2 projection and is infrared finite; and Hq 8 
arises via a twist-4 V2 x twist-2 Vl projection. 

The basic building blocks for annihilation are matrix elements of the operators 
(qib) V - A {q2q3)v-A, {qifyv-Ai&q^v+Ai and {qib) S -p{q2q?>)s+Pi denoted A l { f {h) (ViV 2 ) , 
A^ f{ - h \v x V 2 ), and A% f w (VM) , respectively, see e.g., Figure 3. The first quark bi- 
linear corresponds to the B meson, the superscript % (/) indicates a gluon attached to 
the initial (final) state quarks in the weak vertex, and by convention V 2 (Vi) contains 
a quark (antiquark) from the weak vertex. A3' and A^~ dominate the AS = 1 QCD 
penguin annihilation amplitudes. The latter are expressed as (ViV 2 h \T E h \B)QCD = 
-ic(V 1 V 2 )f B fvJvA{ViV2), where 



ih Cf 



C 3 A\"+ (C 3 + N c C 1 )A f 1 "+ C,Af+ (C 5 + N c C e )A(- h } . (16) 



The arguments (V1V2) have been suppressed. The c coefficients are again determined by 
the vector meson flavor structures. For the electroweak penguin annihilation amplitude 
63 EW , substitute £3,4,5,6 — > £9,10,7,8, respectively. A3' and A^~ arise at twist-3. A%~ is 
given by 

At = -na s /^^(?^0lH-)^ 2 W4-+^^^ 2 W<H-)4-)- (17) 
Jo v m bfv 1 v u m b f v . 2 u l vJ 

For A3' substitute fy^ <-> 4>± — > 4>\\, $6,a — ► and change the sign of the second 
term. The integrals over u and v are logarithmically divergent, corresponding to the 
soft gluon limit u,v — ► 0. For simplicity, the asymptotic distribution amplitudes are 
used, as in [2TJ|22|; non-asymptotic SU(3)p violating effects will be discussed shortly. 
The logarithmic divergences are again replaced with complex parameters, Xa = (1 + 
gAe tipA )\ia ms/Ah, yielding 

— 18 (^# + ^) (2 ^- 3)(1 -^- (18) 

For A f 3 '° substitute f£> «-> and (2X A - 3)(1 - X A ) -> (2A^ - 1)(2 - A^). The 
contributions of A^' - and A3' to the helicity amplitudes are formally of 0(l/m 2 ) (fs 
scales like m b 1 ' 2 ). However, note that as p~ A and p° A are varied from to 1, A%~ and 
A3' increase by more than an order of magnitude. 

A summary of power counting at next-to-leading order is given in |25[ l2"Tj . As ex- 
pected, each quark 'helicity-flip' costs 1/m in association with either one unit of twist, or 
form factor suppression. A ±1 change in vector meson helicity due to a collinear gluon in 
a higher Fock state also costs one unit of twist, or 1/m. In addition, annihilaton graphs 
receive an overall 1/m suppression. An apparent exception is provided by the (twist-3) 2 
contributions to A^; they contain a linear infrared divergence which would break the 
power counting. {A l { 2 would be promoted to (9 (1/m 2 ) but would remain numerically 
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small, as can been by parametrizing the divergence as (m^/A/J^e^ with, e.g., 



K 



3). 



However, the divergence should be canceled by twist-4xtwist-2 effects, see below. Re- 
gardless, (jlj) remains formally true in QCD factorization. The first relation in (|3J) has 
also been confirmed recently in SCET [28] . We expect that the power counting obtained 
in QCD factorization will be reproduced for all corresponding graphs in SCET. 

Amplitudes involving twist-4 vector meson projections remain to be explicitly evalu- 
ated S5J. Twist-4 xtwist-2 projections give rise to Hg 8 . However, these effects should be 
similar in magnitude to (twist-3) 2 contributions to the positive helicity hard spectator 
amplitudes, which were found to be small. The twist-4 xtwist-2 contributions to 
must cancel the non-vanishing (twist-3) 2 contributions, since must vanish by equa- 
tions of motion. This condition leads to new Wandura-Wilczek type relations between 
the products of twist-4 xtwist-2 and (twist-3) 2 light-cone distribution amplitudes. These 
relations should insure cancelation of the aforementioned linear divergence in A l {~ 2 by 
twist-4 x twist-2 effects |25|. Finally, twist-4 x twist-3 projections give rise to A^ + and 
A(' + ; however, these amplitudes should be both formally and numerically suppressed 
by 0{l/m 2 ) compared to A\~ and A3 , respectively. We have also not explicitly con- 
sidered graphs in which higher twist two-body vector meson projections are replaced 
with higher Fock-state projections of same twist containing collinear gluons, e.g., qqg. 
The latter are expected to receive additional suppression at each twist, e.g., 20% [29J. 
These corrections, especially the tree-level twist-3 contributions to the coefficients a^'~ , 
should be included j2H]- However, they will not alter our conclusions, given the large 
uncertainties that have already been assigned to the power corrections. 

Expressions for a few B — > VV amplitudes are given below, 



Ah 



(*F ST^ \(s) Ah „h j_ „v,h 1 h ^„P,h , oh 1 oh 



V2 



p=u,c 



EW 



A h 



p=u,c 



_ \ * \(s) Ah nP> h „V)h \ X oh I Qh i oh 

I 1= A p ^-pK* a 4 n"4,EW "i" °puP2 i* P3 +P3,EW 



Ah 

^B°->p+p- 



,G 



—% 



V2 ^ p pp 



p=u,c 



Ah 

^■B-^p-p 



.Gi 



E ^< 



p=u,c 



h h\ h ^ h 

Opu\^i + a 2 ) + -a3 ; EW + o a 4,EW 



3 
— ( 

2 



(19) 



where P?(ViV 2 ) = b h i {V x V 2 )f B fv x fvjA\ riV2 . For B~ K*-<j> add the term 5 pu f3% to 
Ag ^ K , 0( ^. The arguments (V1V2) of a\ , a\ , ffi have been suppressed, but it is under- 
stood that they are to be identified with the subscripts (V1V2) of the prefactors Ay iV2 . 
The new annihilation coefficients bf, and amplitudes for other decays are given in |25| . 
(To first approximation, all annihilation coefficients except 63 can be ignored in the 
above amplitudes.) b\ 2 arise from current-current operator annihilation graphs. b\ arise 
from QCD annihilation graphs with different flavor topology than 63, and b\ EW are the 
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Figure 3: Quark helicities for H§ 7 (left), A3 (right). The upward lines form V2. 
analogous electroweak annihilation coefficients. 

3 Numerical analysis 

In our numerical analysis the inputs of [22] are used, except for the form factors for 
which we take C|f* (p) = A± .16 (.31 ± .13) and <;f* w = .3 ± .12 (.22 ± .09). The 
expanded ranges reflect the lower values recently obtained for from the lattice (30] 
and from B — > K*^ jHB E2]- We also take r 1 ^ ' p = .05 ± .05, spanning existing model 
determinations |E] |H]- In the evaluation of the hard-scattering and annihilation graphs, 
a s and the Wilson coefficients are evaluated at an intermediated scale ph = y/KfJl, 
with p e [7715/2, m&]. The renormalization scales are varied independently in the three 
classes of graphs: hard-scattering (/xjf), annihilation (p^), and those containing form 
factors. The quantities X4 and Xh parametrizing the logarithmic divergences are varied 
independently for unrelated convolution integrals, with qa,h G [0, 1] and <£a,h G [0, 27r]. 
The default values are g — <p — 0. 

The form factor dependence can be approximately factored out of the (CP-averaged) 
B — > pp branching ratios, yielding 

10 6 (Br p - p o, Br p+p -) = (18.1+|;f, 27.2 j£«) x 

The error bars include the uncertainties due to variation of all remaining inputs, added 
in quadrature. The measured values are 10 6 Br ex ? = 26.4 ± 6.4 (average of [TBI I15j). 
and 10 6 Br^+ p p _ = 30 ± 4 ± 5 |Ej. We also obtain f L (p~P°) = -97tm ± -01 and 
fi(p~P + ) = -^-Ao-m- The first and second sets of error bars are due to form fac- 
tor and remaining hadronic uncertainties, respectively, again added in quadrature. The 
form factor errors are reduced substantially when the branching ratios are constrained to 
lie in their measured ranges, to approximately t'xa f° r P P° an d -'x,l f° r P + P~ |2Sj- The 
relatively small polarization uncertainties are due to the absence of (in p~p°), or CKM 
suppression of (in p + p~) the QCD penguin amplitudes. The measured polarizations are 
/f P (p~P°) = -961$ (average of P5||E|) and fT P (p-p + ) = .99 ± -03±;g| [H]. Thus, the 
predicted longitudinal polarizations are in good agreement with experiment, and with 
naive power counting. 

Averaging the Belle and BaBar B — ► <pK*° (CP-averaged) measurements fHJ Ell 
□H yields = .49 ± .06 and 10 6 Br exp = 10.61 ± 1.21, or 10 6 Br^ xp = 5.18 ± .86 



\V ub \ g 
.0037.31 



(20) 



S 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Pa. Pa 



Figure 4: Br L ((fiK* ) vs. p° A (left), Br T (0i^*°) vs. p A (right). Black lines: default inputs. 
Blue bands: input variation uncertainties. Yellow bands: additional uncertainties from 
logarithmically divergent power corrections. Thick line: Br™ ax , see text. 

and 10 6 Br^ p = 5.43 ± .88, where Br^ = Br_|_ + Bry is the total transverse branching 
ratio. Without annihilation we obtain 10 6 Br L = and 10 6 Br T = .6±;|±:|, 

where the second (first) set of error bars is due to variations of Xh (all other inputs). 
In Figure 4, Br^ and Br^ are plotted versus the annihilation parameters p° A and p A 
entering A3' and A3 , respectively. The black curves are obtained for central values 
of all inputs, with default values for all annihilation and hard spectator interaction 
parameters other than p A , p A . The blue bands are obtained by adding the uncertainties 
due to variations of the inputs in quadrature, keeping default annihilation and hard 
spectator interaction parameters. The widths of the bands are dominated by the form 
factor uncertainties. The yellow bands also include, in quadrature, the uncertainties 
due to variations of all Qh,a-,<£h,a and p^. The thick curve gives the maximum values 
obtained for Br-p under simultaneous variation of all inputs. The absolute branching 
ratios suffer from large theoretical uncertainties, as is usually the case. Nevertheless, 
it is clear that the contributions of A{'° and to the QCD annihilation amplitudes 
can be 0(1) numerically even though they are formally 0(l/m 2 ). This can be traced 
to the quadratic dependence on the divergences (X A ) and the large coefficient N c Cq 
in 63 h . The quantities X A and X A , as well as the renormalization scales and form 
factors entering A and A~ are, a-priori, unrelated. Figure 4 therefore implies that the 
measurements of Br^ and Br t can easily be accounted for simultaneously. According to 
Figure 4, fi((j)K*~) ~ 50% can also be accounted for given that the <j)K*~ and <j)K*° 
amplitudes only differ by a small current-current operator annihilation graph. 

A test for right-handed currents 

In Figure 5 (left) the predicted ranges for f±/f\\ and Br T are studied simultaneaously 
for B -> <fiK* in the Standard Model. The 'default' curve is again obtained by varying 
g A in the range [0,1], keeping all other inputs at their default values, and the error 
bands are obtained by adding uncertainties in quadrature as in Figure 4. Evidently, the 
second relation in (jlj) holds at next-to-leading order, particularly at larger values of Br-r 
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Figure 5: f±/f\\ vs. Biy in the SM (left), and with new RH or LH currents (right). 
Black lines, blue bands, and yellow bands are as in Figure 4. Thick lines: (/±//||) max in 
the Standard Model for indicated ranges of r±. 



where QCD annihilation dominates Br^ and Bry. We also plot the maximum values 
attained for f±/f\\ under simultaneous variation of all inputs. The result is sensitive to 
r_i_, as it largely determines the relative signs and magnitudes of the 'form factor' terms 
in A~~ and A + , see (j2J). The thick black curve (corresponding to Br™ ax in Figure 4) and 
blue curve give maxima for r±_ > 0, in accord with existing model determinations, and 
rj_ > —.10, respectively. A ratio in excess of the Standard Model range, e.g., f±/f\\ > 1.5 
if rj_ > 0, would signal the presence of new right-handed currents. We mention that non- 
vanishing CP-violating triple products in pure penguin decays like B — > <fiK* would not 
be a signal for right-handed currents if significant strong phase differences 0mod7r) 
existed between ^4 ,]| an d A± [HBHH11- There is some experimental indication for such 
phase differences [TU] . which is to be expected if annihilation amplitudes are important. 

Right-handed currents are conventionally associated with effective operators Qi, ob- 
tained from the Standard Model operators Qi by interchanging V — A <-> V + A. The 
final states in ^4 ,|| (-4-l) are parity-even (parity-odd), so that the i'th pair of Wilson 
coefficients enters as [35 j 

Ao,\\cxCf M +C^-Cf p , A ± <xCf u +C? p +C? p . (21) 

The different combinations allow for large modifications to JJ. J'l suffers from pro- 
hibitively large theoretical uncertainties. However, f±/f\\ is a much cleaner observable. 
For illustration, new contributions to the QCD penguin operators are considered in Fig- 
ure 5 (right). At the New Physics matching scale M, these can be parametrized as 

(~) (~) (~) M (~) 

C 4 = C 6 = — 3 C 5 = — 3 C 3 = k . For simplicity, we take M pa M w and consider two 
cases: k = —.007 (lower bands) or k = —.007 (upper bands), corresponding to C^(mb) 

or C^(m b ) pa ,18Cfg(m 6 ), and C^(m b ) or C^(m b ) « .25 Cfg(m 6 ), respectively. 
The default curves and error bands are obtained as in the Standard Model case. Clearly, 
moderately sized right-handed currents could increase f±/f\\ well beyond the Standard 
Model range if r± > 0. However, new left-handed currents would have little effect. 
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Dipole operators versus four-quark operators 

The suppression of dipole operator effects in the transverse modes has important 
implications for New Physics searches. For example, in pure penguin decays to CP- 
conjugate final states /, e.g., B — ► <p(K* — > K s tc°), if the transversity basis time- 
dependent CP asymmetry parameters (S/)j_ and {Sf)\\ are consistent with (sin 2(5) j/^Ksi 
and (Sf)o is not, then this would signal new CP violating contributions to the chromo- 
magnetic dipole operators. However, deviations in (Sf)± or (Sf)\\ would signal new CP 
violating four-quark operator contributions. If the triple-products A° and A\ (33] EH] do 
not vanish and vanish, respectively, in pure penguin decays, then this would also signal 
new CP violating contributions to the chromomagnetic dipole operators. (This assumes 
that a significant strong phase difference is measured between A\\ and A±.) However, 

no n- vanishing A\, or non- vanishing transverse direct CP asymmetries would signal the 
intervention of four-quark operators. The above would help to discriminate between 
different explanations for an anomalous time-dependent CP asymmetry in B — > <pK s , 
i.e., S^Ks, which fall broadly into two categories: radiatively generated dipole operators, 
e.g., supersymmetric loops; or tree-level four-quark operators, e.g., flavor changing (lep- 
tophobic) Z' exchange jSH], -R-parity violating couplings [37], or color-octet exchange 
[3*5] . Finally, a large f±/f\\ would be a signal for right-handed four-quark operators. 

SU(3)f violation and B — ► pK* 

We have seen that the large transverse (j)K* polarization can be accounted for in 
the Standard Model via the QCD penguin annihilation graphs. Would this necessarily 
imply large transverse pK* polarizations? To answer this question we need to address 
SU(3) p flavor symmetry breaking in annihilation. For simplicity, we have thus far esti- 
mated the annihilation amplitudes using asymptotic light-cone distribution amplitudes 
|2U[ • However, for light mesons containing a single strange quark, non-asymptotic 
effects should shift the weighting of the distribution amplitudes towards larger strange 
quark momenta. SU (3) p violation in processes involving ss popping versus light quark 
popping, e.g., annihilation, can therefore be much larger than the canonical 20%, due to 
the appearance of inverse moments of the distribution amplitudes j^H] • This can account 
for the order of magnitude hierarchy between the B — > D°ti and B — > D^K~ rates [39J. 
Similar considerations may also explain the 0(1) flavor violation empirically observed in 
high energy e + e~ fragmentation, e.g., in kaon versus pion multiplicities, or K* versus p 
multiplicities at the Z. In particular, the relative probability for ss popping versus uu 
or dd popping in JETSET fragmentation Monte Carlo's must be tuned to ~ .3 [40J. 

The dominant B — > VV QCD annihilation amplitudes A% ~ and A^'° involve products 
of inverse moments, see (fT7|). The SU(3)f violation discussed above can be estimated by 
including the second and third terms in the Gegenbauer expansions for the distribution 
amplitudes. The first Gegenbauer moments cxi t ±, a\ ii determine the asymmetries of the 
corresponding leading-twist distribution amplitudes, i.e., the inverse moments of 4>± are 
given by (u~ 1 )±, (u~ 1 )± = 3(1 ± a\ t ± + Q2,±), and similarly for <f>\\. Note that the first 
moments vanish for the symmetric <fi and p mesons. For illustration, two sets of intervals 
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Figure 6: Scatter plots of \A~ K ^ {(j)K*) / A~ K ^ {pK*)\ vs. \A~ K ^{<PK*)\ for <p A ^ 0, 
and Ax^(<j>K*)/A- K .fc{pK*) vs. A K ^(<pK*) for <p A = 0, for full (green), halved 
(blue) and default (solid line) ranges of Gegenbauer moments, see text. 

for the moments are considered: aflni) = .2 ± .2, &2L{\\) = -1 ± .3, a 2_i_(||) = ± .3, 
as in [22]; and a more restrictive set with same central values but halved intervals. 
a i ±(||) > is required, since the s-quark should carry the larger fraction of the K* 
light-cone momentum. Similarly, we require that the light-cone momentum fraction of 
the s-quark (light antiquark) in the K* is greater than (less than) that of the quark 
(antiquark) in the p and (ft or, equivalently, 

(OjS) > ^'tV < u_1 >5i) < ^K) > ( 22 ) 

which imposes the constraints ^fl^) > l a fj-(||) — a 2±(||)l- ^he logarithmic divergences 
in the inverse moments are parametrized as before. For simplicity, the Xa are taken 
equal and independent of the final state, with qa < 1 and tpA £ [0,27r]. In j5T], SU(3)p 
violation was studied with asymptotic distribution amplitudes by varying the Xa- 

The scatter plots in Figure 6 illustrate SU(3)p violation in the QCD penguin an- 
nihilation amplitudes A^^P^^R*) and A~ K * ^(pK*), see ()19|) . For simplicity, only 

the contributions of A^~ are included (to good approximation, the other terms in 6j 
can be neglected). Two cases are shown: arbitrary strong phase tfA £ [0,27r]; and 
vanishing strong phase ipA = 0. The Gegenbauer moments are sampled in the inter- 
vals given above, subject to the constraint (J2"2"|) : pfi lies in the usual range, pa € [0, 1], 
and the remaining inputs are set to their default values. For comparison, the default 
non-annihilation (j)K*° amplitude (in units of —iXcGp/V^) is ~ .026, with negligible 
strong phase. The total negative helicity B° — ► <pK* amplitude observed is about a 
factor of three larger, corresponding to Aj.^^^ (K*(f>) ~ .05 and, according to Figure 6, 
.9 < A- K ^{pK*)/A- K ^{(t)K*) < 2.5. We therefore expect / L (iT*°p ± ) ~ fd<pK*°) in 
the Standard Model. The other pK* modes, containing 'tree-level' amplitudes, will be 
discussed in [25]. BaBar has measured f L (p°K* + ) = OMto^t-oxti CHI- Given the large 
errors, this is still consistent with the low cf)K* longitudinal polarization. 
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4 Conclusion 



We have presented an analysis of polarization in B decays to light vector meson pairs 
beyond naive factorization, using QCD factorization. Formally, the longitudinal polar- 
ization satisfies 1 — Jl — 0(l/m 2 ), as in naive factorization. However, we saw that the 
contributions of a particular QCD penguin annihilation graph which is formally 0(l/m 2 ) 
can be 0(1) numerically in longitudinal and negative helicity AS — 1 B decays. Con- 
sequently, the observation of fL(4>K* ~) ~ 50% can be accounted for in the Standard 
Model, with large theoretical errors. However, /l(p + p°) and /l(p + P~) are predicted to 
be close to 1 with small theoretical errors, in agreement with observation. We have shown 
that the ratio of transverse rates in the transversity basis satisfies T±/T\\ = 1 + 0(l/m), 
in agreement with naive power counting. A ratio in excess of the predicted Standard 
Model range would signal the presence of new right-handed currents in dimension-6 four- 
quark operators. The maximum ratio attainable in the Standard Model is sensitive to the 
B — > V form factor combination (1 +my/me)A 1 — (1 —myjm^V or r±, see (J3J), which 
controls helicity suppression in form factor transitions. Existing model determinations 
give a positive sign for r±, which would imply T±(<j)K*)/r\\(<j)K*) < 1.5 in the Standard 
Model. However, the maximum would increase for negative values. The magnitude and 
especially the sign of r± is an important issue which needs to be clarified with dedicated 
lattice studies. 

The contributions of the b —>■ sg dipole operators to the transverse modes were 
found to be highly suppressed, due to angular momentum conservation. Comparison 
of CP-violation involving the longitudinal modes with CP-violation only involving the 
transverse modes, in pure penguin AS = 1 decays, could therefore distinguish between 
new contributions to the dipole and four-quark operators. More broadly, this could 
distinguish between scenarios in which New Physics effects are loop induced and scenarios 
in which they are tree-level induced, as it is difficult to obtain 0(1) CP-violating effects 
from dimension-6 operators beyond tree-level. 

We have seen that the asymmetry of the meson light-cone distributions gener- 
ically leads to 0(1) SU(3)f flavor symmetry violation in annihilation amplitudes, as 
pointed out in j^HJ. In particular, ss popping can be substantially suppressed relative 
to light quark popping. This implies that the longitudinal polarizations should satisfy 
fa^K* ) < f L {(j)K*) in the Standard Model. Consequently, fo^K* ) « 1 would 
indicate that [/-spin violating New Physics entering mainly in the b — > sss channel is 
at least partially responsible for the small values of fL{<f>K*). One possibility would be 
right-handed vector currents; they could interfere constructively (destructively) in the 
perpendicular (longitudinal and parallel) transversity amplitudes. Alternatively, a par- 
ity symmetric realization would only affect, and increase the perpendicular amplitude 
[33] . Either case would lead to T± > T\\, and could thus be ruled out. A more exotic 
possibility is tensor currents; they would contribute to the longitudinal and transverse 
amplitudes at subleading and leading power, respectively. If left-handed, i.e., of the form 
s<7 M „(l + j 5 )bsa^ v {l ± 7 5 )s, then T ± « T\\ would be maintained. Finally, 0(1) SU(3) F 
violation is possible in all QCD penguin amplitudes, given that the annihilation topology 
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components can be comparable to, or greater than the penguin topology components. 
This is especially true of decays to VV and VP final states which, unlike decays to PP 
final states, do not receive large contributions from (S — P)(S + P) chirality penguin 
topology matrix elements. Certain applications of SU(3)f symmetry in B decays should 
therefore be reexamined. 

Acknowledgments: This work originated during the 2003 Summer Workshop on Fla- 
vor Physics at the Aspen Center for Physics. I would like to thank Martin Beneke, 
Alakhabba Datta, Keith Ellis, Andrei Gritsan, Yuval Grossman, Dan Pirjol, and Ian 
Stewart for useful conversations. I am especially grateful to Matthias Neubert for dis- 
cussions of many points addressed in this paper. This work was supported by the De- 
partment of Energy under Grant DE-FG02-84ER40153. 

References 



[i] 

[2] 
[3 

[4] 
[5] 
[6] 
[7] 



[9 
[10 

[11 
[12 

[13 



G. Altarelli et al., Nucl. Phys. B 208, 365 (1982). 

I. Bigi. M. Shifman, N. Uraltsev and A. Vainshtein, Phys. Lett. B 328, 431 (1994). 

J. Charles, A. Le Yaouanc, L. Oliver, O. Pene and J. C. Raynal, Phys. Rev. D 60, 
014001 (1999) |arXiv:hep-ph/9812358| . 

M. Beneke and T. Feldmann, Nucl. Phys. B 592, 3 (2001) |hep-ph/0008255 . 
A. Ali, G. Kramer and C.-D. Lu, Phys. Rev. D 58, 094009 (1998). 



P. Ball and V. M. Braun, Phys. Rev. D 58, 094016 (1998) |hep^p h/9805422] . 

L. Del Debbio, J. M. F lynn, L. Lellouch and J. Nieves, Phys. Lett. B 416, 392 
(1998) |hep-lat/ 9708008 . 

P. Colangelo, F. De Fazio, P. Santorelli and E. Scrimieri, Phys. Rev. D 53, 3672 
(1996) [Erratum-ibid. D 57, 3186 (1998)] |arXiv:hep-ph/9510403| . 

M. Bauer, B. Stech and M. Wirbel, Z. Phys. C 34, 103 (1987). 

B. O. Lange and M. Neubert, |arXiv:hep-ph/0311345 



R. J. Hill, T. Becher, S. J. Lee and M. Neubert, |arXiv:hep-ph/0 404217 



G. Burdman and G. Hiller, Phys. Rev. D 63, 113008 (2001) |arXiv:hep-ph/0011266| . 

M. Beneke and T . Feldmann, Nucl. Phys. B 685, 249 (2004) 
|arXiv:hep-ph/0311335|. 



[14] I. Dunietz, H. R. Quinn, A. Snyder, W. Toki and H. J. Lipkin, Phys. Rev. D 43, 
2193 (1991) 



14 



[15 
[16 
[17 
[18 
[19 

[20 

[21 

[22 
[23 
[24 

[25 
[26 

[27 

[28 
[29 

[30 

[31] M. Beneke, T. Feldm ann and D. Seidel, Nucl. Phys. B 612, 25 (2001) 
|arXiv:hep-ph/0106067| . 

[32] A. Ali and A. Y. Parkhomenko, Eur. Phys. J. C 23, 89 (2002) 
|arXiv:hep-ph/0105302| . 

15 



J. Zhang et al [BELLE Collaboration], Phys. Rev. Lett. 91, 221801 (2003) 
|arXiv:hep-ex/0306007| . 

B. Aubert et al. [BAB AR Collaboration], Phys. Rev. Lett. 91, 171802 (2003) 
|arXiv:hep-ex/0307026| , and |hep-ex/0303020| 

B. Aubert et al. [BABAR Collaboration] , |arXiv:hep-e x/0404029 , submitted to Phys. 
Rev. Lett. 

K.-F. Chen et al. [BELLE Collaboration], Phys. Rev. Lett. 91, 201801 (2003) 
|arXiv:hep-ex/0307014| . 

J. G. Smith [BABAR Collaboration], talk at Rencontres de Moriond QCD04, March 
2004, |arXiv:hep-ex /04060 63| A. Gritsan [BABAR Collaboration], L BNL seminar, 
April 2004 |http://costard.lbl.gov~gritsan/RPM/BABAR-CQLL-0028.pdf| . 

M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Phys. Rev. Lett. 83, 
1914 (1999) hep-ph/9905312| ; Nucl. Phys. B 591, 313 (2000) |hep-ph/0006l"24| ; 

M. Ben eke, G. Buchalla, M. N eubert and C. T. Sachrajda, Nucl. Phys. B 606, 245 
(2001) |arXiv:hep-ph/0104110| . 

M. Beneke and M. Neubert, Nucl. Phys. B 675, 333 (2003) |arXiv:hep-ph/0308039| . 

H. Y. Cheng and K. C. Yang, Phys. Lett. B 511, 40 (2001) |arXiv:hep-ph/0104090| . 

X. Q. Li, G. r. Lu and Y. D. Yang, Phys. Rev. D 68, 114015 (2003) 
larXiv:hep-ph/0309136| . 

A. L. Kagan, in preparation. 

P. Ball, V. M. B raun, Y. Koike and K. Tanaka, Nucl. Phys. B 529, 323 (1998) 
lhep-ph/9802299| . 



A. L. Kagan, talks at Super B Factory Workshop, Honolulu, January 2004; Les 
Rencontres de Physique de la Valle d'Aoste, February 2004. 

C. W. Bauer, D. Pirjol, I. Z. Rothstein and I. W. Stewart, |arXiv:hep-ph/0401188| 

A. Hardmeier, E. Lunghi, D. Pirjol and D. Wyler, Nucl. Phys. B 682, 150 (2004) 
|arXiv:hep-ph7030717l| . 



D. Becirevic, arXiv:hep-ph/0211340 



[33] G. Valencia, Phys. Rev. D 39, 3339 (1989). 
[34] A. Datta and D. London, arXiv:hep-ph/0303159 

[35] A. L. Kagan, talks given at SLAC Summer Institute, August 2002; SLAC 10 36 
Workshops, May and October 2003. 

[36] V. Barger, C. W. Chiang, P. Langacker and H. S. Lee, Phys. Lett. B 580, 186 (2004) 
|arXiv:hep-ph/0310073| . 

[37] A. Datta, Phys. Rev. D 66, 071702 (2002) |arXiv:hep-ph/0208016|; B. Dutta , 
C. S. Kim and S. Oh, Phys. Rev. Lett. 90, 011801 (2003) |arXiv:hep^ph/0208226|. 



[38] G. Burdman, |arXiv:hep-ph/0310144| 



[39] S. Mantry, D. Pirjo l and I. W. Stewart, Phys. Rev. D 68, 114009 (2003) 
|arXiv:hep-ph7 0306254| . 

[40] T. Sjostrand, |arXiv:hep-ph/9508391| 



[41] M. Beneke, eConf C0304052 (2003) F0001 |arXiv:hep-ph/0308040|. 



16 



